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The microscopic structure of an isolated vortex line in superfluid neutron matter is studied by 
solving the Bogoliubov-de Gennes equations. Our calculation, which is the starting point for a 
microscopic calculation of pinning forces in neutron stars, shows that the size of the vortex core 
varies differently with density, and is in general smaller than assumed in some earlier calculations 
of vortex pinning in neutron star crusts. The implications of this result are discussed. 
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The interior of a neutron star constitutes the only 
known physical system close to infinite nuclear and neu¬ 
tron matter. Most efforts to describe this system have 
mainly focussed on the role of nucleonic interactions at 
zero and hnite temperature, while less attention has been 
paid to a microscopic description of excited states of the 
system, including superfluid vortex lines induced by the 
rotational state of the star. On the other hand, there 
are important observables of astrophysical relevance that 
might be influenced by the presence of vortex lines. A 
mutual check of the nuclear many body physics and of 
the theory of neutron star interiors comes for example 
from the study of pulsar glitches, sudden increases in the 
spinning frequency of the crust of the pulsar, followed 
by a slower tendency to conditions close to the original 
ones. It is thought that glitching events represent a di¬ 
rect manifestation of the presence of superfluid vortices 
in the interior of the star ||], the triggering event being 
an unbalance between the hydrodynamical forces acting 
on the vortex and the force of interaction of the vor¬ 
tex with the nuclei present in the crust (pinning force). 
There have been quite large uncertainties regarding the 
value of the pinning force, leaving room for quite oppo¬ 
site views about the validity of the vortex pinning model 
l^. One source of uncertainty is related to the value 
of the pairing energy gap in uniform neutron matter, a 
quantity strongly dependent on the value of the neutron 
particle-particle matrix element near the Fermi surface 
1^-^. A second problem is due to the very approximate 
way of treating vortex states in neutron matter. Usu¬ 
ally a vortex is seen as a cylinder of normal matter (the 
vortex core) of radius equal to the BCS coherence length 
^0 ~ O.Sh^kF/‘2mA where kp is the Fermi wavenumber, 
m is the neutron mass and A is the neutron pairing gap 
l^-H. The pinning of the vortex to the nucleus, also 
treated as a classical object, is due to the loss of super¬ 
fluidity that occurs when the two objects superimpose. 

Several experiments and theoretical calculations in 
type II superconductors show that a vortex is a quite 
complicated quantum object |^. Bound states that can 
be formed in the center of the core as revealed by the in¬ 
crease in the density of states (investigated by scanning¬ 
tunneling microscopy) may change the local pairing gap 
in the core and the vortex size. In addition, they show 


that superfluidity is not completely suppressed in the 
core. The latter result can be in principle predicted 
also with the Ginzburg-Landau theory (GL), but unfor¬ 
tunately the GL theory is not reliable far from the transi¬ 
tion temperature, which is the case occurring in neutron 
stars. 


In the present Letter we propose to study the structure 
of a vortex in superfluid neutron matter using a micro¬ 
scopic, fully quantum-mechanical approach. More specif¬ 
ically, we shall make use of the Bogoliubov-de Gennes 
equations 0 , that have been successfully employed 
to study vortices in type II superconductors 00 
and more in general non-homogeneous superconductivity. 
The use of a model well-tested in the laboratory repre¬ 
sents, in our view, a major merit of this kind of approach. 
There are however important differences between vortices 
in neutron matter and those in superconductors as stud¬ 
ied in Ref. that make a direct application of the 

previous findings impossible, and require a new ab-initio 
calculation. The first obvious difference is in the basic 
constants setting the scale of the problem, like the mass 
of the particles and the density regime. The difference 
is not merely quantitative, since in our case the param¬ 
eter ^okp defining the importance of quantum effects is 
smaller than in the type-II superconductors, resulting in 
the possibility of strong deviations from previous solid- 
state calculations and also from semiclassical methods. 
Secondly, due to the charge neutrality of neutrons, mag¬ 
netic fields do not play any role in our calculation and 
vortices are generated by the rotational state of the star. 
A further difference is that the temperature in supercon¬ 
ductors is an important variable of the problem, while in 
our case, due to the high value of the Fermi temperature 
compared to the interior temperature of a neutron star, 
it can be set equal to zero. A fourth major difference 
is the role of the inter-particle interaction, which in the 
case of neutron matter is strongly density-dependent. Fi¬ 
nally, we shall be mostly interested in information about 
the pinning force, which depends in a delicate manner on 
the length scale of the vortex core. 


A vortex state in neutron matter can be described by 
the Bogoliubov-de Gennes (BdG) equations 0: 
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+ W{r) - Ef^ U^{v) + A(r)F,(r) = E,U^{v) (1) 

-(^^ + W{v) - Ef^U,{v) + ^*{v)V^{v) = E,V,{v) (2) 

where p = —ihV, U and V are the quasiparticle am¬ 
plitudes, Ef is the Fermi energy, A(r) is the (space- 
dependent) pairing gap, m is the neutron mass, W{v) 
is an external potential due to e.g. a nucleus, and the 
subscript i represents all relevant quantum numbers. We 
write the quasiparticle states as 

Uiir) = j^Un^ip) exp[i{p - 1/2)0] exp[i{kzz)] (3) 

Vi{r) = exp[i{p + l/2)9]exp[i{k^z)]. (4) 

Here L is the length of the cylinder, taken equal to 1 in 
our calculations since all quantities are per unit length 
of the vortex, n is a radial quantum number and p is 
half an odd integer. The quantity kz is the (conserved) 
momentum in the z direction. The pairing gap has to be 
calculated self-consistently as 

A(r) = -5 ^ C/.(r)b/*(r) (5) 

i;0<Ei<hQ 

where the sum is over quasiparticle states having energy 
Ei smaller than a cutoff and g > 0 is the pairing 
strength. The meaning of these parameters is quite dif¬ 
ferent for neutron matter compared with the solid state 
case, and will be discussed later on. Due to the above an¬ 
gular dependence of the quasiparticle states, the pairing 
gap has the form A(r) = A(p) exp[—10], corresponding 
to a vortex with one quantum of circulation. Follow¬ 
ing 1^] we expand the quasiparticle states in terms of 
cylindrically symmetric Bessel functions and impose the 
boundary condition of zero gap at the edges of a cylin¬ 
der of radius R that is chosen sufficiently large. More 
specifically, the basis functions are chosen as 

<Pjm{p) = yyy , / = 1; • ■ ■ j A^j (6) 

where m = /r ± ^ is an integer and ajm is the jth zero 
of the Bessel function The dimension N of the 

basis is chosen large enough to ensure convergence and 
stability of the quantities of interest. The quasiparticle 
amplitudes are expanded as 


(P) =^Cnj4’j^,-^ip) 

j 

(7) 

(p) 

(8) 
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For a given value of p Eqs. (|^J|) can then be written as 
a 2N X 2N matrix eigenvalue problem 


where the superscript T denotes the transpose of a ma¬ 
trix, (ct^i ,..., Cyjjv, djii ,..., djiFf ), 


T± 


2m* 




Ekl-k 


2 

F 


Sii 


and the matrix A is given by 
R 

4‘3p.-l/2{p)Mp)4‘3'F+iMp)pdP- 

We have in these equations set the mean field W (r) equal 
to zero, corresponding to an isolated vortex line in infi¬ 
nite neutron matter. Interactions between neutrons set 
up an internal mean field which can be simulated by re¬ 
placing the neutron mass m by an effective mass m*. 
However, for dilute neutron matter m* « m |Q|. Start¬ 
ing from an appropriate approximation to A(p), we solve 
Eq. (^ for several values of p to obtain the quasiparticle 
amplitudes U and V. A new approximation to A(p) is 
then obtained from Eq. (||). The procedure is iterated 
until convergence is achieved. 

The BdG equations as written above have been spe¬ 
cialized to the case of a point-like neutron-neutron inter¬ 
action of the form n(r, r') = —gS{r — r'). This is the same 
approximation as used in Refs. llQ-h^ and should be con¬ 
sidered even more reliable in the case of neutron matter, 
where the range of the inter-particle interaction is of the 
order ~ 1 fm, and is thus larger than the length scale of 
the inhomogeneities, which is of the order of several fm. 
In dealing with a zero-range force, a cutoff energy hfl 
has to be introduced for the gap equation to converge. 
The density dependence of the pairing gap as calculated 
with finite-range forces, see for example [^|, can be easily 
translated to a dependence of g on the density. We shall 
thus fix hid = 50 MeV and use the value of g that repro¬ 
duces the gap at large distance from the vortex core as 
calculated for a homogeneous system. 

Fig. 0 shows the pairing gap as a function of the dis¬ 
tance from the vortex axis, A(p) using the above formal¬ 
ism and for three different values of the neutron Fermi 
wavenumber. In all cases the gap increases from a value 
zero to an asymptotic one Aoo- The latter quantity rep¬ 
resents evidently the value of the gap as calculated for 
uniform neutron matter. In the figure we have fixed the 
value of the effective pairing strength g to reproduce the 
value as calculated from the bare Bonn A potential [||. 
Like for vortices in type H superconductors, the gap is 
found to increase linearly for small values of p before 
reaching the asymptotic value. The linear rate of increase 
of the gap defines a coherence length that for consistency 
with ref. ||^ we call ^ 2 , given by 



f T- A 
A^ T+ 


4 '„ = En'i’n 


lini A(p) = Aoo-p. 
p^o 42 


(9) 
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( 10 ) 








As visible from the figure this coherence length repre¬ 
sents an appropriate length scale of the vortex core. A 
second length scale, the BCS coherence length can be 
defined as the distance from the vortex axis at which the 
energy gain of the Cooper pair due to the velocity field 
becomes of the order of the energy necessary to break 
the pair (or also by the condition that the superfluid 
velocity Vg = Ti/mp coincides with the Landau velocity 
for the suppression of superfluidity, 2/S.oa'm/Tikp)^ and it 
is found ^0 = 16.88fci7’/Aoo. This BCS coherence length, 
which microscopically represents the size of a Cooper pair 
[p^ , is the one usually employed to calculate the interac¬ 
tion of a vortex with the nuclei in the crust of a neutron 
star &I0- In particular, matter is supposed to be 
completely normal (A(p) = 0) for p < Co and superfluid 
(A(p) = Aoo) for p > Co- In a pictorial view, this length 
can be obtained intersecting the gap at infinity, Aoo with 
the curve g{p) = 0.814 x fi^kp/’^mp. The resulting gap 
profile is shown in Fig. || for the cases kp = 0.1 fm~^ 
and kp = 0.5 fm~^. It is evident that this definition 
can be a reliable estimate of the vortex core size only 
when the intersection occurs for values of p smaller than 
C 2 . The plot for kp = 0.1 fm“^ is a case where the in¬ 
tersection occurs where the gap has already saturated. 
According to the common view of neutron star pinning, 
all the matter at p < Co should be considered in a nor¬ 
mal state. This is not the case, as clearly visible from the 
figure. For kp — 0.5 fm“^. Fig. || shows that Co < C 2 , 
and the approximation above is probably not too bad, al¬ 
though it somewhat overestimates the gap in the region 
Co < P < C2- 

Table | shows Co and C 2 as functions of the Fermi mo¬ 
mentum, where the pairing strength g has been fitted to 
the ^Sq gap of the Bonn A potential Q. It is visible from 
the table that the value of C2 seems to be uncorrelated 
with Co- la general, C 2 remains within values between 6 
to 10 fm, while the values of Co are much more scattered. 
Due to the fact that Friedel oscillations may slightly dis¬ 
tort the vortex structure at small distances from the core 
p^ , we also studied some vortex length scales defined 
without making use of the behavior close to the axis, for 
example the distance from the vortex core at which the 
gap is a fraction of Aoo, typically 50 % and 90 %. As 
expected, these quantities are larger than C2 , but show a 
similar trend. Overall, we found these lengths to agree 
better with C 2 than with Co- 

To understand the scaling of the vortex size as a func¬ 
tion of the Fermi momentum, we also performed calcu¬ 
lations keeping the pairing strength fixed and changing 
only the Fermi wavenumber, even if this produced un¬ 
realistic gaps. We found C2 to be a decreasing function 
of the wavenumber, a result which agrees qualitatively 
with previous work on vortices in superconductors [OJ^ 
where the vortex size is predicted to scale as ~ , but 

is very different from the behavior of the BCS coherence 
length, which scales proportionally to kp/A^a- In the 
realistic case, where the pairing strength is fitted to the 
gap in uniform matter, the coherence length ^2 increases 


as a function of Thus, the density dependence of the 
pairing strength g is also very important, and the vor¬ 
tex size appears to depend in a quite complex way both 
on g and on the Fermi wavenumber. Studying the dis¬ 
tribution of the eigenvalues, we have found indeed that 
several bound states are formed in the vortex core. This 
feature is different from the solid-state case 111 where 
essentially only one bound state at each value of the an¬ 
gular momentum p is formed. The value of ^2 is prob¬ 
ably affected in a complex way by all the bound states, 
making a simple scaling law difficult to work out. As a 
partial conclusion, however, we can state that when vary¬ 
ing the Fermi wavenumber and the pairing strength over 
the range of values relevant for ^Sq neutron pairing, the 
coherence length ^2 does not change by more than 50 %, 
while ^0 can vary by a factor of ten. 

We are now able to speculate on the applications of 
our results to the case of vortex pinning in the crust. 
The use of ^2 instead of ^0 as the length scale of the 
vortex core, implies that the pinning regime is approx¬ 
imately the same in the whole inner crust. In particu¬ 
lar the super-weak pinning, proposed in the case of very 
large vortex core and in which several nuclei may be en¬ 
veloped on one single plane of the vortex, is at variance 
with our findings. A second consequence already pre¬ 
dictable on the basis of the data presented is that the 
pinning force exerted by a nucleus on a vortex should 
scale like ^ A^kp/^ 2 - We expect an increase of the pin¬ 
ning force of the order ^ C 0 /C 2 compared with the cases 
where ^0 is used for the size of the vortex core, and this 
can be a very large factor for low densities. 

We also examined the distortion of a vortex line in 
the presence of a mean field potential of cylindrical sym¬ 
metry. To simulate the presence of a nuclear cluster, 
we chose a Woods-Saxon shape for the potential W{r) 
with nuclear parameters of the order of those expected 
in the inner crust of a neutron star. In presence of 
a potential field, the diagonal parts of the Hamilto¬ 
nian matrix become ^ -|- where the ma¬ 

trix elements of the potential matrix are of the form 
lT,y =/dpp(/>j^±i/2(p)lT(p)((>y^±i/2(p)(5jy. We found 
a decrease of the coherence length in presence of the field. 
This shows that the shape of the vortex line is sensitive 
to the presence of the nucleus, an effect that might have 
important consequences in vortex pinning. This in fact 
contrasts with the assumption implicit in many calcula¬ 
tions of vortex pinning in neutron star crusts, where the 
vortex core is assumed to be inert during the pinning 
state. 

A further advantage the BdG formalism, not exploited 
in the present study, is that the Hamiltonian is diagonal¬ 
ized, resulting in a consistent calculation of the energy, 
while in other approaches the local density approxima¬ 
tion jl^ or Ginzburg-Landau approaches need to be 
used. The calculation of the pinning energy and force 
of a nucleus with a superfluid vortex, technically more 
involved but conceptually straightforward, is thus under 
consideration Ip3. 
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kp (fm 

Aoo 

^0 (fm) 

6 (fm) 

0.1 

0.06 

26.80 

11.46 

0.2 

0.38 

8.73 

6.59 

0.5 

1.91 

4.42 

10.20 


TABLE I. Gap at infinity and the length scales and ^2 as functions of density. 
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FIG. 1. Pairing potential A(p) as a function of p for an isolated vortex in neutron matter. 
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FIG. 2. The exact pairing potentials at kp = 0.1 and 0.5 fm ^ compared with those obtained by using the ^0 cutoff 
prescription described in the text. 
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